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Abstract 

We reconsider the problem of the solid-liquid-vapour contact-line on a disordered substrate, 
in the collective pinning regime. We go beyond scaling arguments and perform an analytic 
computation, through the replica variational method, of the fluctuations of the line. We show how 
gravity effects must be included for a proper quantitative comparison with available experimental 
data of the wetting of liquid helium on a caesium substrate. The theoretical result is in good 
agreement with experimental findings for this case. 

1 Introduction 

When a liquid partially wets a solid, the liquid-vapour interface terminates on the solid, at the 
contact line. If the solid surface is smooth, then at equilibrium, we expect no distortions of the 
contact line, and the Young's relation giving the contact angle in terms of the interfacial tensions 
holds, that is 

7si> - 7s/ = 7cos(6'eg) (1) 

where 'y = "fiv, and 9eq is the equilibrium mean contact angle. 

We consider a case where the substrate is weakly heterogeneous and where the heterogeneities are 
"wettable" defects, leading to a space dependance of the interfacial tensions jsv and 7^;. Favoured 
configurations are those where the liquid can spread on a maximum number of defects. We thus 
expect distortions of the contact line which tends to be pinned by the defects. Moreover, the energy 
due to the liquid-vapour interface induces an elastic energy of the line. The competition between 
the elastic energy and the pinning due to the disorder gives rise to a non trivial wandering of the 
line, a typical example of the general problem of manifolds in random media |^ . The case of the 
contact line is of special interest for several reasons. There exists by now good experimental data 
for the correlations which characterize the wandering of the line Q. On the theoretical side, the 
problem presents two specific features. The elasticity of the line is non local. The pinning energy 
due to the surface heterogeneities is, up to a constant, a sum of local energy contributions due to 
the wetted defects. It has therefore non-local correlations which are of the "random field" type in 
the usual nomenclature of manifolds in random media. 
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In this paper we will consider the case of collective pinning where the strength of the individual 
pinning sites is small, but pinning occurs due to a collective effect. This seems to be the relevant 
situation for the experiments. The case of strong pinning by individual impurities was studied by 
Joanny and De Gennes ||5|. Collective pinning is a particularly interesting phenomenon since the 
balance between the elastic energy and the pinning one results in the existence of a special length 
scale first discussed by Larkin in the context of vortex lines in superconductors Q. This Larkin 
length is such that the lateral wandering of a line, thermalized at low temperatures, on length scales 
smaller than ^, are less than the correlation length A of the disorder (range of the impurities), while 
beyond ^ the lateral fluctuations become larger than A and the line probes different impurities. 
The Larkin length scale diverges in the limit where the strength of disorder goes to zero. At zero 
temperature, the line has a single equilibrium position when its length is smaller than ^, while 
metastable states appear only for lengths larger than ^. Therefore one can think of the contact 
line, qualitatively, as an object which is rigid on small length scales (less than ^), and fluctuates on 
larger length scales. A third length scale which is relevant for the discussion is the capillary length 
Lc, which is the length scale beyond which effects due to gravity become important: the line then 
becomes "flat" in the sense that its fluctuations do not grow any longer with the distance. 

The collective pinning of the contact line was first addressed by Vannimenus and Pomeau [0. 
They considered the case of very weak disorder in which the Larkin length ^ is larger than the 
capillary length. So their analysis only probes the "Larkin regime" of length scales less than ^, in 
which there exist only very few metastable states. A more complete qualitative picture, making 
clear the role of ^, can be obtained by some scaling arguments originally developped for some related 
problems by Larkin P] and Imry-Ma For the case of the contact line, these arguments were 
introduced by Huse [|lO| and developped by De Gennes Q and by Joanny and Robbins 1^]. They 
lead to interesting predictions concerning the lateral fluctuations of the line: these should grow like 
the distance to the power 1/2 on length scales less than ^, and to the power 1/3 on larger distances 
on length scales between ^ and Lc- More recently, Kardar and Ertaz [11| have performed a dynamic 
renormalisation group calculation for the contact-line at zero temperature, subject to a uniform 
pulling force and also find a roughness exponent 1/3. These scaling laws have been confirmed in 
recent experiments on the wetting of helium on a caesium substrate Q, confirming the validity of 
the collective pinning picture in this case. 

The aim of our paper is to go beyond the scaling analysis and provide a quantitative computation 
of the correlation function of the line. We use the replica method together with a gaussian variational 
approximation, with replica symmetry breaking []l2| . This approach, which is exact in the limit 
of large dimensions, is known to give good results even for one dimensional systems as this one 



14, 15]. It confirms the scaling exponents derived before, but also provides the prefactors and 



a full description of the crossover between the two regimes around the Larkin length. 

The paper is organised as follows: We introduce the model in section (^). In section (|^), we 
present for completeness a scaling argument which gives the roughness exponents, and we obtain an 
expression for the Larkin length by a perturbative approach. In section (Q), we present the replica 
calculation and compute within a variational approximation the full correlation function in the limit 
of low temperatures, first neglecting gravity effects, and then including them. In section (|5|), we 
compare our theoretical prediction with experimental data. 
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Figure 1: Sketch of the experimental set up 



2 The model 



Consider a situation given by figure (|l|), where the hquid wets an impure substrate which is sUghtly 
inchned with respect to the horizontal. We denote by (x, y) the space co-ordinates of the substrate. 
The excess energy per unit area due to pinning is given by 



e{x,y) = -isi{x,y) - ^svix.y) - lsi{x,y) - 7sv{x,y) 
resulting in a total pinning energy 



(2) 



dx / dy e{x,y) (3) 
/o Jo 

where is the height of the the contact line at the abscissa x (overhangs are neglected), and L 
the width of the substrate. As for the pinning energy per unit area or force per unit length e{x, y), 
we shall suppose that it is gaussian distributed, which is the case if it results from a large number of 
microscopic interactions, and that it has local correlations on length scales of order A. Specifically, 
we choose 



W 



e{x,y)e{x',y') = -^S{x 



x')C 



y-y 



A 



(4) 



where the correlation function C(r) is normalised to C(0) = 1 and C"(0) = — 1, and decreases fast 
enough to zero for r ^ 1. The asymmetry introduced in (^) between the two directions x and y 
is for computational convenience. In most physical situations, the distribution of disorder should be 
isotropic in the x — y plane, leading to a correlation in the x direction on length scales of order A. 
As we shall explain below, we have found that this correlation has only small effects, which exist 
only on very short distances and are not relevant experimentally. As for the shape of the function 
C(r), we shall first use for simplicity 



C{r) = f{r) = exp(-rV2) , 



(5) 



and we shall later comment on the modification of our result for more general correlations. 

We must also add to the random potential term, a capillary energy term, which, if we neglect 
gravity and suppose that the slope of the liquid-vapour interface varies smoothly, is given by 



cap-^ / , ^\k\\m\' 

2 ./^<|fc|<^ 2Tr 



(6) 
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with c = ^ — Oeq being the average equihbrium contact angleQ. 

The final hamiltonian is thus 



H = ^ [ ^\k\ + r dx Vix, ^x)) 

2 J^<|fe|<^ ivr Jo 



(7) 



where V{x, <!>)=/ dy e{x, y). As a sum of independent gaussian variables, V{x, $) is gaussian 
distributed, and up to a uniform arbitrary random shift we can choose: 

V{x,^)V{x',<^') = -W5{x-x)f(^( ^~^ ^ ^ (8) 



where f{u) is a function which grows as for large |n| . Its precise form depends on the correlation 
function C of the energy per unit area, and is given in the simple case (||) by 

f{u^) = \u\ / dv e-"'/2 _ (1 _ g^«V2) (9) 



This model provides a good description of the problem of a contact line on a disordered substrate 
under the following hypotheses: 

• The slope of the liquid-vapour interface is everywhere small. 

• The length of the contact line is small compared with the capillary length, so that one can 
neglect gravity. 

• The defects in the substrate are weak and give rise to collective pinning. 

The main part of our work will be dedicated to the analytical study of this simplified model. We 
shall then examine the corrections due to gravity, to more general correlations of the disorder, and 
to the correlations in the x direction. The quantity which is measured experimentally and which 
we shall compute is the correlation function of the position of the line: 



mx)-^{y)Y) (10) 

where thermal averages are denoted by angular brackets and the average over disorder by an overbar. 
As we shall see, in different length regimes, this correlation increases as a power law, which defines 
locally the wandering exponent C, from: 



I2? . (11) 

3 Perturbation theory and scaling arguments 

For completeness we rederive in this section an expression for the Larkin length by perturbation 
theory, and review the scaling derivation of the roughness exponents. 

3.1 The Larkin length 

On a sufficiently small length scale, we can assume that the difference in heights between any two 
points is small compared with the correlation length A of the potential. We can thus linearize the 
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potential term such that V{x, $(x)) ~ V{x, 0) — e{x)^{x). This leads to a random force problem 

W 

with a force correlation function e(x)e(x') = -ri^Six — x'). Rewriting the hamiltonian as: 



c f dk 



we get for T — > and A <^ \x — x'\ <^ L, 



m - 4# 



^ 1 r dk |e(A;)|2 
~ Yc 2^ \k\ 



(12) 



— -— ^, 2W f dk (1 - cos(k(x - x'))) W , 

The wandering exponent in the Larkin regime is given by C = 1/2- The hnear approximation is 
no longer valid when |$(x) — <l>(x')| becomes of the order A. Typically |x — x'| is then of order 
c^A^ 

^ = ^ , where ^ is the so called Larkin length. 

3.2 Roughness exponent for large fluctuations 

On length scales larger than the fluctuations of the line are greater than the correlation length 
A and perturbation theory breaks down. One can estimate the wandering exponent by a simple 
scaling argument as follows ||l2| . The hamiltonian is given by (0) and we can no longer linearize the 
potential term in (^). 

We consider the scale transformation, x — > Ix, ^>(x) l'^^{x), V{x,^{x)) — > l'^V{x,^{x)). 
Imposing that the two terms in the hamitonian scale in the same way and that the potential term 
keeps the same statistics after rescaling, we have 



A = 2C - 1 and 2A = -1 + C (14) 

and so C = 1/3. Note that this is less than the value 1/2 obtained in the Larkin regime. On a still 
larger length scale (larger than the capillary length), we expect the line to be flat and C = 0- 

This exponent can be recovered by the following Imry-Ma argument j^, |l^, 0, |9[. On a scale 
L, the line fluctuates over a distance ^. The elastic energy contribution then scales as c<&^. As for 

, — /l$ 

the pinning energy, since it is a sum of independant gaussian variables, it scales as VWAy 

where y/WA is a measure of the pinning energy on an area A^ and — ^ an order of magnitude of 

A^ 

the number of such pinning sites. Minimising the total energy — VWA^ with respect to 
we get <I> ~ A ( — I with ^ ~ 



U W 



4 The replica computation 

4.1 Computation of the free energy 

We now turn to a microscopic computation of the free energy F = —TlnZ. Since the free energy 
is a self averaging quantity, the typical free energy is equal to the average of F over the disorder. 
We compute it from the replica method with an analytic continuation of Z", for n ^ [^]. The 
n^^ power of the partition function 
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a=l { a a •'^ ) 



gives after averaging over the disorder 



a=l 



(16) 



where 



a,b 



(17) 



We note that the expression of the free energy is invariant with respect to a translation of the 
centre of mass of the Hne ^CM = x dx ^{x) = j^^i^ ~ 0)- centre of mass so that 

there is no integration on the A: = mode. The partition function cannot be computed directly. 
Fohowing [^], we perform a variational calculation based on the variational hamiltonian 



Ho 



dk 
2^ 



(18) 



a,b=l 



where G ^ is a hierarchical Parisi matrix. 
The variational free energy 



pn n 



gives up to a constant term, 



— = lim — 

L n^o n 



where 



and 



du 



(19) 



(20) 



(21) 



1 r dk 

Bab = ^j^ {Gaa{k) + G,,{k) - 2Gab{k)) 

The optimal free energy is obtained for a matrix G verifying the stationarity conditions 



(22) 



which read 



dGab 



0, 
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Gab = -^^f ) fo'^ « ^ ^ 

(23) 



More details on this approach can be found in ||, [l^]. 
4.2 The replica symmetry breaking solution 



To solve equations (pSj), we suppose that the matrix G has a hierarchical replica symmetry breaking 
stucture d la Parisi. We can write G~i^ = {c\k\ —d')5ab — c^ab- is thus parametrised by a diagonal 
part c\k\ — fj, and a function (y{u) defined on the interval [0, 1] 0. The optimisation equations for G 
can then be written as: 

^M = ^/\^J (24) 

with 

o [-Ah 

Bi^) = ^ J ^m)-9ik,u)) (25) 
The solution to these equations is described in appendix A. It is best written in terms of the function 



a]{u) = ua{u) — I dv a{v) (26) 



which is given by 

(27) 

W 

H('") = — XT for u>Uc 

We shall give the value of the breakpoint u^, in the (experimentally relevant) limit of low tempera- 
tures. Defining 

Tc = vrc— , (28) 

we get Uc — T /T(.. From the expression 

du [c](u) 



1 + 



n2 [a][u)+c\k\ 



(29) 



we obtain in the regime A ^ |x — x'| <^ L 

^< ($(x) - $(x'))2 > = Ah(^^^^ (30) 
^From equations (Eq), we know that the off-diagonal elements of G^^ do not depend on k. 
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where 



and ^ 



4 /x\ 



2/3 



3 V vr j 



dk- 



(l-cos(A:)) /-(^M)'/' dw 2x 



+ 



u)^ + 1 Svr 7o k{x/'K + k) 



(31) 



The function Ti is the analytical prediction in the regime where gravity effects can be neglected. It 
is plotted in figure ^ and has the following asymptotic behaviour. For small x, 7i{x) ~ \/\x\ and 
for large x, TC(x) ~ 1.14|x|^/^. Therefore the predictions in the various scaling regimes, including 
pref actors, is as follows. 
When |x — x'l <C ^, (Larkin regime): 



< (<I>(x) - <i>(x'))2 > ~ A 
When |x — x'l ^ ^, (random manifold regime): 



x — X 



1/2 



< (<I>(x) - $(x'))^ > ^ 1.14A 



X — X 



1/3 



(32) 



(33) 



Before turning to the comparison with the experiment, we first compute various correction factors 
to this formula. 



4.3 A more general form of the disorder 

We show that even in the more general case where we only impose that the correlation function 
of the potential has the asymptotic behaviour /(|m|) ~ \/M" for large \u\, the height correlation 
function can be put in the form 

^< mx) - cI>(x"OF^ = Ag(^^^^ (34) 

in the limit T — > 0, and for A <C |x — x'| ^ L. The derivation of Q, for an arbitrary correlation 
function /, is given in appendix B. 

4.4 Effect of the cut-off A at small x 

Coarse-graining the force in the x direction on scales of order A leads to a discretized (in x) version 
of (Q), which is equivalent to the form which we use, but with a cut off of order A at small x. On 
scales comparable to A, there are thus corrections to equation (|3l|) due to this cutoff, which are 
easily computed. When we take into account the cut-off A, the correlation function TC becomes 



(35) 

where ^ ~ '^^^ effect of the cut-off is to shift slightly downwards the theoretical curve, especially 
in the region of small x. 
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0.0 I ' ' ' ' 1 

0.0 5.0 10.0 15.0 20.0 25.0 

X 

Figure 2: The theoretical predictions for the rescaled height correlations. From top to bottom, the 
curve 7i{x) (|30|) where gravity is neglected (A = tt^/Lc = 0), and the curves Tlg{x,X) for values of 
A such that A = 0.1, 0.2, 0.4, 0.6, 0.8. Notice the need to have a very small ratio of Larkin length to 
capillary length in order for the zero gravity result to be good. 

4.5 Effect of gravity 



In the geometry considered, the effective capillary length is given by ./ , where a is the tilt 

y pg sin a 

angle of the substrate with respect to to h orizontal Q |. To take into account gravity we must replace 
the kernel |A:| in the hamiltonian by ^yk^-\-l/L^, where Lc is the capillary length. Generalising 
the previous calculation we have computed the corrections due to gravity. The fluctuation of the 
line at distance x now depends on two parameters, x/^ and A = tt^/Lc- The computation done 
in appendix C expresses [(t](m) in the terms of an inverse function Z~^. When the Larkin length 
is sufficiently small compared with the capillary length, that is when ^ is less than approximately 
0.5Lc, we have in the limit where T <C Tc 
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[o"](u) =0 for u < ui 

[a]{u) = fic I^^ ( ( — ) ) for ui < u < u'^ 
\n\UcJ J 

/r /,/\3/2x (36) 

with ui = Uc[ — 2r(0) I and ~ — 

T 

As for u'^, it is slightly larger than and also of order — . When the capillary length goes to 
infinity ui tends towards and u'^ towards Uc- The function I{x) is given by 

= ( r \ 2 / dk \ -^/^ 

\Jo {^/^+i + x)y \ Jo i^/^Ti + xy) 

and for large x, T{x) ~ x. This asymptotic behaviour ensures that we do recover the results of 
section (|5|) when goes to zero. The correlation is now given by 



X — x' vr^ 



< {^{x) - ^{x')f > = AW, ( ^— , ) (38) 



where 



4 p (l-cos(A;x/7r) r^<l''-^"^ dw \l~^{w^/\) 



i Jo Vk^ + X^ J(. 



2uc f°° „ (1 - cos(A;x/^) XI-^(l/X) 
+-— / dk 



(39) 



nl{x,X) c,x\- dw —I-\w'/\) + -^XI-\l/X) (40) 



3<yo + A2 AJ-i(l/A) + VA:2 + A2 

The asymptotic behaviour of is different from that of 7i. For small x 

'{AI{0))i/3 

and for large x,7ig{x, A) tends towards a constant depending on A. 

Prom the previous equations we can see that gravity has a significant effect when A = — becomes 

of order 1, where ^ is the Larkin length and Lc the capillary length. Moreover we can also note 
that the correction for small A to the case without gravity is of order A^/^. The limit A going to 
is thus a rather slow one. This is illustrated in figure (pi). 



5 Comparison with experiment 
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5.1 The experimental set up 



We have fitted the data from experiments carried out by C.Guthmann and E.Rolley Q with our 
theoretical curve. The experiments study the wetting properties of Uquid hehum 4 on caesium 
below the wetting transition temperature which is about 2K. Above that temperature caesium is 
wetted by helium. In the experiments carried out by Guthmann and Rolley the substrate consists 
of caesium deposited on a gold mirror which is slightly inclined with respect to the horizontal (see 
figure ^). The wetted defects are small areas on the substrate where the caesium has been oxydised. 
The experiments are carried out on a range of temperatures going from about IK to 2K. There 
is a constant inflow of helium at the bottom of the helium reservoir to maintain the contact angle 
to its maximum value 6a, the advancing angle, which is in general different from the equilibrium 
contact angle 6(,q (see figure [l|). This is necessary because otherwise, the liquid would recede and 
the contact angle would shrink to zero due to strong hysteresis. Height correlations are calculated 
from snapshots of the advancing line when it is pinned. The incoming helium is regulated to ensure 
that the line moves with a small velocity and so we can probably suppose can we are just at the 
limit of depinning. 

The predicted order of magnitude of Tc given by (|2|) is ^ = ^7^^^^ {Oe,)^\ g.^^ 
the impurities can be measured experimentally and is of the order of 20fj,m. We thus expect the 
correlation length to be a few times this size. Its precise value depends on the details of the disorder. 
For temperatures not too close to transition temperature, 9eq ~ 20 degrees, and 7 ~ 10^^ Kin~^ . 
This leads to a typical estimate Tc ~ 10^ K, in the experimental conditions of 0]. Therefore T/Tc 
is of order 10~^ and the system is effectively at low temperatures, justifying the low temperature 
limit in our computations. 



5.2 Comparison between experimental data and the theory neglecting gravity 

We consider experimental data for height correlations Tiexpi^) — \l < {^{x) — <I>(0))^ > for temper- 
atures T = 1.72, 1.8, 1.9, 1.93-ftr. The equilibrium angle 6eq, the liquid-vapour interfacial tension 7 
and the "potential stength" W depend on the temperature, and so the various experimental curves 
correspond a-priori to different values of the Larkin length. Instead, the correlation length A which 
depends only on the substrate, is expected to remain constant. We have thus fitted the experimental 
curves to the theoretical prediction in the absence of gravity (^) and (^8|), with the same A but 
different ^'s. 

We proceed as follows. We first minimise the error function 

1 ^ V 

^ = E E ^n^i/^j) - ^e.pixi) (41) 

2^ Nj V / 

j 

with respect to A and for j G [1,4], where j denotes a given experimental curve at the temperature 
Tj, Xi the experimental points, Nj the number of points of curve j, and the correlation length at 
temperature Tj. This procedure yields A ~ 18/Lim and £ ~ O.Sfim?. 

In figure (^), we show the four experimental curves, together with the corresponding theoretical 
fit ATC{x/^j). The values of the Larkin lengths are given in table 1 below. The dependance of 
on the temperature is related to the variations of W, 7 and 9(.q, which are not known well enough for 
a detailed comparison between theory and experiment. It should be noted that the determinations 



11 



40.0 - 



.ra 20.0 - 




0.0 



0.0 



500.0 1000.0 
Transverse distance in microns 



1500.0 



Figure 3: The experimental data for the fluctuation of the hue as a function of the distance (circles). 
From top to bottom, the temperatures are 1.93K, 1.9K, 1.8K and 1.72K. The full curves are the 
theoretical predictions in absence of gravity (|30|) , using the values of Larkin lengths shown in table 
1. 



T(K) 


1.72 


1.8 


1.9 


1.93 




220 


145 


92 


53 




4.1 


2.7 


1.7 


1 



Table 1: Values of the Larkin length scale used in the flt to the experiment shown in flgures (^Q). 



of A on the one hand, and the on the other hand, are strongly correlated. (To give an idea, 
with A = 19/im, = 60/im, ^2 = 104/im, ^3 = 163//m, ^4 = 248//m the error (41) differs from the 
previous case by about only 10%). In the absence of detailed information on the experimental error, 
it is thus difficult to give an error bar on ^. On the other hand, the ratios of the Larkin lengths in 
different experiments are much less sensitive to this correlation. They are given in table 1 as well. 

In figure (^, we show a collapse plot of all the experimental curve on the theoretical one with 
gaussian correlated disorder. We have rescaled each experimental curve by the corresponding 1/^ 
in the x direction and by 1/A in the y direction. This collapse gives several interesting results: the 



curves nearly collapse one onto the other, as expected from the general form (34). Furthermore 
it seems that the simplest correlation function that we have studied in most details (|5|) gives a 
reasonable fit to the data. We have also checked that the small x cut-off discussed in section ( [4.4[ ) 
is irrelevant. However there is also clearly, a systematic difference at large distances, which we shall 
now discuss. 
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Figure 4: Collapse plot of all the experimental data on fluctuations of the line at various tem- 
peratures onto the theoretical prediction neglecting gravity (|30|), using the values of Larkin lengths 
shown in table 1. 

5.3 Comparison between experimental data and the theory including gravity 

In figure (Q), we note that each of the rescaled experimental curves lies below the theoretical 
curve at large distances. Moreover the experimental curves have a slightly larger curvature than the 
rescaled theoretical curve. These are indications that gravity cannot be totally neglected. Indeed the 
effective capillary length in the experimental conditions is of the order of 2mm, and experimentally 
the correlations are measured for distances up to about 1.5mm, which is actually not small compared 
with the capillary length. 

To check whether gravity has or not a significant effect, we have tried to fit the experiments 



with the full theoretical prediction including gravity (38). Encouraged by the results of the previous 
analysis, we keep to the case of a gaussian correlation function of the disorder given by (§),(|9|). 
We have carried out the same analysis as in the previous subsection, using as theoretical input 
TCg{x, A) instead of Ti.{x). The capillary length Lc is not adjustable: it is calculated for the different 
temperatures from the experimental measurements of 7, and are given in table 2. Therefore this 
new fit has the same number of adjustable parameters as the previous one. We now minimise the 
error function: 

= E E AWg(x,/e„ <,/L^J - He^pix.) (42) 

2^ Nj j=i i=i V / 

j 

As one could expect form the rather slow convergence of the theoretical curves Ti.g{x,\) towards 
the gravity-free one Ti.{x) = Tlg{x,0) at small A (see figure |2|), we find rather different values for the 
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T(K) 


1.72 


1.8 


1.9 


1.93 


Lc inm) 


1855 


1838 


1819 


1823 




480 


425 


365 


295 




1.6 


1.4 


1.2 


1 



Table 2: Values of the experimental capillary lengths, and of the Larkin length scale used in the fit 
to the experiment shown in figure (^). 

parameters. In this case, A ~ 75/um and £g ~ 0.26/xm^. The values of the Larkin lengths for the 
different temperatures are given in table 2. We have the same problem of correlations between the 
determination of A and the Larkin lengths as before. The ratios of the Larkin lengths in different 
experiments are also given in table 2. 



60.0 




0.0 I ' ' ' ' ' ' 1 

0.0 500.0 1000.0 1500.0 

Transverse distance in microns 

Figure 5: The experimental data for the fluctuation of the line as a function of the distance (circles). 
From top to bottom, the temperatures are 1.93K, 1.9i^, 1.8K and 1.72i^. The full curves are the 
theoretical predictions in the presence of gravity (^), using the values of Larkin lengths shown in 
table 2. 

In figure (|5|), we show the experimental curves together with the corresponding theoretical curve 
TCgix, A) including gravity. The fit is clearly better in this case since we have got rid of the systematic 
drift from the theory for large x. The exponents of the two experimentally observed regimes are 
correctly predicted by a theory without gravity, but to obtain the correct values of the correlation 
length and Larkin lengths it is necessary to take into account the effect of gravity. 
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6 Discussion and perspectives 



Our analytic computation for the wandering of a contact line on a disordered substrate fits quite 
well the experimental data. We have shown that gravity effects are far from negligible and should be 
taken into account in order to extract the relevant parameters (which are basically the correlation 
length of the disorder as well as the Larkin length scale) from the experiments. Taking into account 
gravity in the theory clearly improves the fit to the data. Moreover, the theory including gravity 
will allow larger experimental length scales in the analysis. 

A few comments about the validity of our computation are the following. First of all, we have 
used in most of our analysis a simple form for the correlation function of the disorder (P) which is 
not necessarily the correct one. It is true that some of our predictions are independant of this form, 
like for instance the existence of a scaling behaviour in absence of gravity. It is an experimental 
problem to have a better description of the pinning disorder at work, and we have shown that our 
computation can be extended to any type of correlation. The simplest one seems to give already 
a good account of the data. The variational method which has been used in this work is also an 
approximation. So far there is no other analytic quantitative method available, and it seems to 
work quite well, confirming previous evidence found in other problems [13, 0, 15 1. 

A more interesting question concerns our assumption of an equilibrium situation. We have 
supposed that the line is in thermal equilibrium in order to write the usual partition function, and 
at the end of the calculation, we have taken the temperature to zero since, as we can see from the 
numerical values of the parameters, thermal fluctuations are irrelevant. In doing so, we retain only 
the states with the lowest energy. Experimentally the state of the line has no apparent reason to be 
a low temperature equilibrium state. It is difficult so far to characterize fully the metastable states 
that can be reached dynamically by the experimental procedure. They might have generically the 
same statistical properties as the ground state, this is actually under study (fl^). A less interesting 
but may be more realistic alternative, could be that because of gravity, the fluctuations of the line 
do not go much beyond the correlation length. So, even though we are out of the Larkin regime, we 
are not yet deep in the random manifold regime and there are thus not so many metastable states. 

A purely dynamical computation could also be done. While the properties of the unpinned line 
have already been studied in [11|, the dynamics of the pinned object could also be very interesting. 
If the correlation length of the disorder can be made much smaller, so that thermal fluctuation 
are no longer negligeable, we expect the onset of some ageing dynamics [19|. It will be interesting 
both to compute its properties and to measure them. In particular, this system might present a 
nice situation to measure the fluctuation-dissipation ratio in a system which has a full (continuous) 
replica symmetry breaking, offering a chance to measure directly the function [a] (u) . 
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A Computation of the function [cr](w) 

We give in this appendix some details of the calculation of the function [a] (u) . We have from section 

^(^) = ^/{^) (43) 



where 



and from [O], 



9 f rik 

B{n) = -p J ^m)-9ik,u)) (44) 



g{k)-g{k,u) = ———^-——-f ^ , , \ . (45) 
u{c\k\ + [a\{u)) 7„ c\k\ + [a\{v) 



Differentiating (|4^) gives 



^ W = {u)f ( ^ ) (46) 



and replacing B\u) in ( |4^) by its expression 



B'(„) = 

pvrc [ctJ (u) 



leads to 

0"'(u) =0 or 



4^ 1 >/ij(n)\ 
vrcA^ [ct]('u)-^ V 



(48) 



We express ~^2^ in terms of [o"](u) by inverting the second equation of (^8[). This gives: 



Differentiating (|49|), and using expression (46) to express B'{u), and the fact that [o']'(^i) = ucj'{u), 
we get 

r (B{u)\ PttcA^ 3T, 

u = -—u (50) 



where 



/'(x) = , (51) 



From (pop and (|5l|), we have 



18 



vrcuA^ TcU ' ' 

Multiplying both sides of equation ([47|) by u, and using (|52|) we get after integrating over u 

[a]{u)=Au^/^ (53) 
The breakpoint Uc, above which the solutions to ( |5^ and (|5^) are no longer valid is given by 

T T 

When — goes to 0, Uc — — . Now since B{u) tends to infinity when u tends to 0, we have o"(0) = 0. 

J- r. -Lr 



For u > tic, B{u) = B{uc) and [c](n) = [cr](uc). To obtain A, we can differentiate ( |53[ ) and compare 

W 1 

the result with the expression for a'{u). We find A = — — . . For u > Uc, cr'{u) = and so 
W 

[a]{u) = [(j]{uc) = — TT- 
vrcA 

B General form of the correlation function for arbitrary disorder 

In this appendix we derive the height correlation function for a more general form of the function 
/ appearing in the correlation function of the disorder (^). We only impose that /(|w|) ~ -^■\/N" 
for large u where A is some constant, such that f{u) ~ -\/|ti[. We shall keep the same notations as 
in appendix A. In this more general case, equations (p3D,(|47|),(|50|) from appendix A are still valid. 
We define the function for positive 

h'Hx) = Lj^ (55) 

3 

where h~^(x) ~ for large x. The asymptotic behaviour of h(y) for small and negative y is then 

2x 

3 

. We now express B{u) in terms of h, from equations (^) and (|55|). This gives 

B{u)=A^h(-^^u] (56) 



Differentiating the previous equation (|5^) gives 



which can be rewritten as 



2T 



B'{u) = -^h!\-^u] (57) 



2T d 

— -\o^\a\(u) = w}^{-w) (58) 



TTCA^ . . 

where w = u. Integrating (|58|) gives 
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[a]{u) = [o-](e)exp 



3u/2uc 



3e/2uc 



dw wh' {—w) 



For e < ti ^ ?ic, we can use the asymptotic form of h in (5£), which then reads 

/ \ 3/2 

H(^) = H(e) - 



e / 



Now for small u, equation (H^) becomes in this case 



, ^ W (T,\^'^ 1 



and so for small n, since [(7](0) = 0, we get 



7rcA4 

A comparison of this last expression with (|60|) gives 



3/2 



w ( 



3/2 



Replacing this last expression in (|^) and taking e to zero leads to 

W 



\a\{u) 



S{u) 



where 



(59) 



(60) 



(61) 



(62) 



(63) 



(64) 



S{u) = ( ^ 



3/2 



exp 



'iu/2uc 



dw w[ h'{—w) 



2tt;2 



(65) 



For the sake of simplicity, we will suppose that in this case the break-point up to which expression 
( |65|) is valid, is also Uc in the limit of low temperatures. This implicitely requires that /i(— 3/2) = 0. 
Then for u > Uc 



1-3/2 / 3 \ 

S{u) =S{uc) = exp J dw w\h'{-w) - 



(66) 



When /' has the simple form (51), h i—w) = — ^, and we recover the expression of [<7l(n) derived 

2W''' 

in appendix A. In the limit of low temperatures and for A ^ |x — x'| <C L, the height correlation 
function is given by 



< ($(x) - ^{x')f > = AG 



X — X 



(67) 



with 
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X 4/x\2/^ /-^ „(l-cos(A;)) /-(^M)'^' 



1 



f + exp 



dtfl h'{—w) 



2w' 



2^ /"^^r 



(1 - cos(A;)) 



k^x/n + kexp^— J dw ^h' {—w) — -^-^ 



(68) 



C Effect of gravity 

In this appendix, we consider a specific case of the disorder given by equations (§) and (P). To take 

into account gravity, we must replace the kernel \k\ by ^/\k\^~+Ji^ , with fi = —, where Lc is the 

Lc 

capillary length. The equations derived in appendix A are thus no longer valid. If cr'{u) is not zero 
then, equation ( |49| ) of appendix A becomes 

where 



IC{x)=cf^g{^] (70) 



with 



Differentiating (|4^) gives 



dk 



g{x) Jo (VFTT + x)2 



(/C'(H(u)))3/2 



(71) 



2a^(^) 1 
vrpc /C([crJ(uj) 

Differentiating (|69|), and using dH) and (H), we have 

1 I B{u) _ T 1 

A2 reu/C(H(n)) ^ ^ 

Differentiating the previous expression ([73|), and using (|7^), we can express [o"](n) as 

/C(H(U)) _^^3/, ^^^^ 



where A is a constant to be determined. Replacing in equation (46), B'{u) by its expression ( |69D 
and reexpressing (^Tj) using ( [7^ ) leads to 



21 



^ W 1 {ulC'{[a]{u))fl^ 



Comparing the previous expression with ( |7^ gives 



W 1 

To express [cr] , it is convenient to introduce the function X 

which is strictly positive and increasing. X(0) ~ 0.87 and Tix) goes as x when x goes to infinity. 
The inverse function is thus defined on the interval [X(0),oo]. Since must be continuous 

and [(t](0) = 0, the function \a\ necessarily has a first plateau where [i7](n) = from n = up to a 
value u\ given by 

For ui < u < ti^, where u'^ is the new break-point to be determined 

H(„) = c,i-(|(!l)''^) (79) 



The break-point u'^ is obtained using ( |7^ ) and is given 



B(u') 1 2 f dk 



A2 ^</C'(H«)) 3 

and when T goes to zero 



n,=.</C'(H«)) =<5'(M^) (81) 



where 



(82) 



(V^^ + n)3 \Jo (Vg2Tl + n)2 
Since (7' is a strictly increasing function, with ^'(O) = S/ir'^ and g'{oo) = 1, in the limit of low temper 

atures Uc < u' < / . Since I(x) is almost linear, we can suppose that [(7](n' ) ~ — —7 ( — ) 
g [0) TrcA*\UcJ 
and so equation (|8T|) can be rewritten as 



^ = ^'^'g'm (83) 



with 
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n=U'^X^^ and A = ^ (84) 

\\UcJ Lc 

We can solve for u'^ perturbatively using the solution without gravity. As a first approximation, we 

can take for A the its value in the absence of gravity. We can then solve numerically for and for 
u'^. For our experimental data, we get by this method u'^. ~ Uc- For u > u'^, [cr]{u) = The 
height correlation function is then given by 



where 



' ' 3 Jo VA;2 + A2 7(AX(o))i/3 >X-\w^/\) + VP + A2 

3^ ^ Jo VFTA2 AX-i(l/A) + VFTA2 



(85) 



(86) 
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